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ABSTRACT 

We present Chandra observations of the galaxy cluster A4059. We find strong evidence that the FR-I radio 
galaxy PKS 2354-35 at the center of A4059 is inflating cavities with radii ~ 20 kpc in the intracluster medium 
(ICM), similar to the situation seen in Perseus A and Hydra A. We also find evidence for interaction between the 
ICM and PKS 2354-35 on small scales in the very center of the cluster. Arguments are presented suggesting that 
this radio galaxy has faded significantly in radio power (possibly from an FR-II state) over the past 10 s yr. 

Subject headings: galaxies:jets, galaxies:clusters:individual (Abell 4059), hydrodynamics, X-rays: galaxies: 
clusters 



1. INTRODUCTION 

Clusters of galaxies are complex dynamical structures and 
their cores are subject to an array of interesting physical pro- 
cesses. Constraints from imaging X-ray observations suggest 
that the hot X-ray emitting intracluster medium (ICM) in the 
core regions of rich clusters is radiatively cooling on timescales 
shorter than th e life of the cluster, giving rise to cooling flows 
(Fabian, 1994, and references therein). The central dominant 



features ar e seen in CXO observations of Hydra A (M cNamara 
et al „ 2000;pavid et al„ 2001|) and Abell 2052 ( Blanton et al., 
2001). 

In this Letter, we present CXO observations of the rich 
galaxy cluster Abell 4059 (z = 0.049). The cD galaxy of 
A4059 hosts the FR-I radio galaxy PKS 2354-35. A short 
ROSAT High Resolution Imager (HRI) observation of this 
source suggested the presence of t wo ICM cavities at the sam e 



galaxy present in many clusters often hosts a radio loud active 
galactic nu cleus (AGN). It has been suggested (e.g., B inney & 
Tabor, 1995) that cooling flows and central cluster radio galax- 
ies are intimately related via complex feedback processes. It is 
easy to see how radio galaxy activity resulting from black hole 
accretion can be associated with a cooling flow. However, the 
impact of a radio galaxy on its environment is much less clear. 

Theoretically, we expect radio jets to inflate cocoons of rel- 
ativistic plasma that expand into the surrounding ICM (e.g. 
Begelman & Cioffi, 1989j|Kaiser & Alexander, 199% R eynolds 



1998 



position angle as the radio lobes ( Huang & Sarazin, 
hereafter HS). In § 2 we discuss our observations, confirm- 
ing the presence of these cavities, and show that A4059 dis- 
plays significant additional morphological complexity. Con- 
straints on models for this source are discussed in § 3, § 4 
presents our conclusions. We assume a Hubble constant of 



et al., 2001b, , hereafter RHB). The energy input by this process 
has recently come un der investigation for i ts potential role in 



heating cluster cores (Briiggen et al., 200 1| ; Quilis et al., 2001 



Reynolds et al., 200k). However, while our simulations sug- 



Hq = 65kms Mpc and qq — 1/2, giving a linear scale 
of 1 kpc arcsec -1 = 0.492 kpc pixel -1 . 

2. observations 

2.1. Data Reduction 

The galaxy cluster Abell 4059 was observed with the Ad- 
vanced CCD Imaging Spectrometer (ACIS) on 24-Sept-2000 
(22.3 ksec exposure) and on 4-Jan-2001 (18.4 ksec exposure). 
The radio nucleus of PKS 2354-35 was placed 1 arcmin from 
the nominal aim point of the back-illuminated S3 chip, placing 
the bulk of the emission on chip-S3. In this Letter, we only use 
data collected on this chip. The data was read out at the standard 
3.2 sec frame rate, telemetered using faint mode, and filtered on 
ASCA event grades. The energy range was restricted to the 0.3- 
10 keV band and corrected for exposure and vignetting, spectral 
fitting was restricted to 0.3-8.0 keV. 

Some of the Sept-2000 data were affected by a period of rel- 
atively high background, with up to 10 counts per second per 
chip for the back-illuminated chips. However, most of the clus- 
ter emission is well above this rate and we decided to use the 
entire data set for imaging. For the spectral analysis, we filtered 
the data on the counts from chip-S 1 to fall within a factor of 1 .2 
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gest that about half of the energy injected by the jets can be 
thermalized in the cluster center, numerical simulations of this 
process still carry a large degree of uncertainty, since limited 
computational resources require significant simplications. In 
order to verify the validity of the assumptions and to design 
future models, we require guidance from observations of radio- 
galaxy/cluster interactions. 

Imaging X-ray observatories, such as the Chandra X-ray Ob- 
servatory (CXO), provide a direct probe of this interaction. 
Both ROSAT and CXO observations of Perseus A have found 
X-ray cavi ties coincident with the radio lobes (B oehringer 
et al., 1993; pabian et al., 2000| ), surrounded by X-ray shells 
which appear to be slightly cooler than the unperturbed ICM 
(see, for example, RHB for a possible explanation). Similar 
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of the quiescent background rate, which rejected 33% of the ob- 
serving window. We then used the quiescent background files 
by Markevitch 7 for background subtraction. The background 
subtracted total count rate on chip-S3 is 5.1 counts/second. 

2.2. Morphological Appearance 

Figure la shows an adaptively smoothed image of the central 
3 arcmin o f A4059. Overlayed is the 8 GHz image by (T aylor 
et al., 1994, , hereafter T94). While the large scale emission ap- 
pears smooth (probably due to the larger smoothing length used 
by the csmooth algorithm), the core of A4059 is not relaxed: 
it is double peaked with one peak at the cluster center and the 
other 15 arcsec south-west of the center. The smoothed image 
reveals a distorted hour-glass like structure (yellow regions in 
Fig. la), centered on the nucleus and oriented perpendicular to 
the radio axis. 

The most interesting detail of the X-ray image are the two X- 
ray holes already noticed by HS, clearly visible in the 5-sigma 
smoothed image. The higher quality CXO image confirms that 
these are real local brightness minima, not just a visual effect 
caused by a bright central bar perpendicular to the cluster elon- 
gation. A noticeable offset exists between the nucleus and the 
axis between the cavity centers, with the nucleus being shifted 
to the south-west by ~ 12" from this axis. The cavity centers 
are roughly 50" apart. 

To assess the significance of the cavities we extracted an az- 
imuthally averaged radial surface brightness profile of the cen- 
tral region, excluding the cavity regions shown in Fig. la (a by- 
eye approximation) and calculated the expected surface bright- 
ness in the two cavities. Based on this estimate, the NW and SE 
cavities are significant to 26 and 5.3 sigma, respectively. How- 
ever, because the central region of A4059 is strongly perturbed 
(roughly inward of 1 arcmin from the center), it is difficult to 
make a rigorous statement concerning the significance of the 
cavities this way. 

However, the outer regions of the cluster appear relaxed. As 
an alternative method, we took radial surface brightness profiles 
from the combined SE- and NW quadrants and the combined 
SW-NE quadrants. As a first order correction for cluster ellip- 
ticity we shifted the radial SW-NE profile outward by a factor 
of 1.13, producing a good match at large radii with the SE-NW 
profile. Because the temperature at large r is relatively uni- 
form, we used an isothermal /3-model to fit the surface bright- 
ness, which represents the outer cluster very well. Because the 
inner cluster regions are not well represented by a /3-model, we 
only used points further than 60 arcsec from the center for the 
fit (r c ~ 50", = 0.52, x /dof = 59/43 = 1.4). 

We then compared the flux measured inside the cavity con- 
tours in Fig. la to that expected from the /3-model. The NW 
cavity has 33 sigma significance compared to the best fit /3- 
model, while the less pronounced SE cavity has 13 sigma sig- 
nificance. Forcing the /3-model normalization to be consistent 
with the observed cavity flux increases the reduced chisquare 
of the fit by a factor of 2.7 (NW cavity) and 1.9 (SE cavity), 
corresponding to a significance of 8.5 and 6 respectively. 

The radio overlay in Fig. la shows that the 8GHz radio lobes 
are only partly coincident with the cavities. The NW lobe cov- 
ers a good fraction of the NW cavity. The SE lobe is much 
smaller than the northern lobe and not spatially coincident with 
the SE cavity. We note, however, that low-frequency radio ob- 
servations sensitive to spatial scales of 0.1 — 1 arcmin, which 



might reveal the full extent of the radio lobes, do not yet exist. 

The central peak contains interesting sub-structure (Fig. lb., 
smoothed to a S/N of 3-4). The brightest sub peak, which is 
well resolved by CXO and has a diameter of 3-4"^ 3-4 kpc, 
is coincident with the core of PKS 2354-35 (within the CXO 
pointing accuracy) and could be emission from the hot inter- 
stellar medium of the central galaxy. The three sub-peaks are 
located aroun d a local brightness minimum. Our hardness ra- 
tio analysis (§2.3) suggests that this minimum is due to lack of 
emission rather than intervening absorption. The close corre- 
spondence between these peaks and the SE radio lobe suggests 
that this substructure might be caused by on-going interaction. 

2.3. Spectroscopic Properties 

We used the adaptive binning routine of Sanders & Fabian 
(2001) to produce hardness ratio maps. We define the hard- 
ness ratios hx = (1 - 2keV)/(0.3 - 1 keV) and h 2 = (2 - 
10 keV)/(l — 2 keV) as the ratio of the counts in the respective 
energy bands and c\ as the surface brightness in the 0.3-10 keV 
band in cts cm -2 s _1 pixel -1 . For ICM observations, h\ is 
mainly sensitive to absorption variations, whereas h 2 is a tem- 
perature diagnostic. While there is no sign of varying absorp- 
tion across the field (as seen in hi), the h 2 map clearly shows a 
radial temperature gradient, as would be expected in a cooling 
flow cluster (Fig. 2). It also shows a global temperature gradi- 
ent, with the S W half of the image appearing hotter than the NE 
half. 

We have extracted a global spectrum of A4059 from the 
central 90 arcsec. Fitting this spectrum with a two tem- 
perature wabs*zwabs* (mekal+mekal ) thermal plasma 
model (Galactic neutral hydrogen column fixed to A^h,g = 
1.45 x 10 20 cm -2 ) results in the best fit parameters kT\ = 
1.34±g;f|keV, kT 2 = 3.90±^|keV, Z = 0.60±g;^, AT H ,z = 

5.36±£;gxl0 20 cm -2 , x 2 / d °.f = 1-29(3 sigma error bars). A 
wabs*zwabs* (mkcflow + mekal) cooling flow model 



provides a similarly reasonable fit (kTi 
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+0.75 
-0.1 



kT 2 = fcT mekal = 3.80t° i^kcV, Z = 0.71 



! +0.09 
-0.09' 



keV, 
M = 

27.6±^ Mq yr-\ iV H . z ="5J54t°-*; x 10 20 cm-\x 2 /dof = 
1.26). 

3. DISCUSSION 

Like Perseus A and Hydra A, A4059 shows clear signs of 
interaction between the radio galaxy and the ICM. However, 
unlike in these systems, there is no one-to-one correspondence 
between the radio lobes and the X-ray holes. We argue be- 
low that the cavities were indeed created by the radio activity 
(given the similarity in position angle and the fact that there are 
no other viable models for producing large ICM cavities), but 
that the radio galaxy has faded since they were formed. This is 
particularly interesting since CXO has recently found clusters 
with cavities but without a radio source (e.g., M cNamara et al., 
2001): A4059 could be a 'missing link' between cavities with 
and without detectable radio lobes, supporting the notion that 
cluster cavities can be created by radio galaxies without having 
to show current radio activity. 

3.1. Evolutionary state 

In the hydrodynamic simulations of RHB, we identified three 
evolutionary phases. Early on, the cocoon is highly overpres- 
sured, driving a strong shock into the ICM. Simulated X-ray 



7 See http://asc. harvard . edu / cal / Link s/Acis/acis/Cal_prods/bkgrnd/ cur rent 
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maps show a cavity surrounded by a thin, hot shell. Once the 
cocoon comes into pressure equilibrium with the ICM, the side- 
ways expansion of the cocoon becomes sub-sonic and the shock 
becomes a compression wave, though active jets can keep the 
advance speed of the jet heads supersonic for much longer, driv- 
ing "sonic booms" into the ICM. For roughly a sound cross- 
ing time, the well defined cavity created during the supersonic 
phase will survive and be observable. Then, hydrodynamic in- 
stabilities will destroy the cavity altogether. 

A strong shock would show up as a sharp feature in our X-ray 
images, which have an effective resolution 8 of 4-5 arcsec at the 
edges of the cavities. The absence of such a shock indicates that 
the source expansion is no longer highly supersonic. We iden- 
tify the hourglass-like feature in the core with the compression 
wave ("sonic boom") found in our hydrodynamic simulations 
(Fig. 1, RHB; in our simulated X-ray maps the brightest emis- 
sion also tends to be in the equatorial plane). 

3.2. Source Power 

We can estimate the radio s ource parameters b ased on the 
presence of the X-ray cavities (|Heinz et al. 1998 RHB, Chu- 
razov et al. 2000). We use the two-temperature fit of § 2 to 
estimate the physical parameters of the ICM. Taking the hot 
emission to arise uniformly in a sphere of 90" radius yields an 
electron density of rihot > 0.009 cm~ 3 . Assuming the cold 
gas is in pressure equilibrium with the hot gas gives an electron 
density of n co id > 0.031 cm~ 3 and a volume filling factor of 
5 x 10~ 3 . 

We assume that both X-ray cavities are completely evacu- 
ated by the lobes and estimate their size from the smoothed 
images by approximating them as spheres. While this is clearly 
a simplification, it will be sufficient for this order of magni- 

8 After adaptive smoothing to a signal to noise of 5. 



tude estimate. A 'by eye' fit of the cavities gives bubble radii 
of rbub ~ 20" (20 kpc). The pressure in the hot phase is 
Picm PS 1-1 x 10 -10 erg cm -3 , relatively close to the minimum 
energy pressure in the lobes of pme ~ 3 — 5 x 10~ 11 ergs cm -3 
(T94). At a minimum, the radio galaxy has to perform "pdV" 
work against the ICM. Including the internal energy of the 
plasma within the cavities, this gives an integrated energy out- 
put of Etot £ 8 x 10 59 ergs. 

The radio galaxy had to inject this energy into the bubbles 
before they floated out of the cluster core. This buoyancy 
timescale is approximately twice the sound crossing time of the 
relevant region of the cluster, r ~ 4?"bub/cs ~ 8 x 10 7 yr, 
where we have used the sound speed for a 4keV gas, c s ~ 
1000 km s" 1 . The time-averaged source power needed to pro- 
duce the cavities is then Lkm ~ E^/t J> 3 x 10 44 erg s _1 . 

Alternatively, we can estimate the source age from the "sonic 
boom" arguments of RHB. The hour-glass structure through 
the cluster center is roughly 50" long (i.e., 25" ~ 25 kpc 
on either side of the center). Following RHB, we equate this 
to the distance traveled by a shock/compression wave which 
moves at least at the sound speed of the hot ICM. This gives a 
source age of r <; 2.4 x 10 7 yrs and a time-averaged power of 

ikin £ E h /r ~ 10 45 ergs" 1 . 

A third, X-ray independent way to estimate the source power 
is based on the radio flux. The flux densities at 5 GHz and 
8 GHz are given b y T94 as 76mJy an d 34mJy, while the 
1 .4 GHz NVSS flux Condon et al. (1998b is 1 .3 Jy. The NVSS 



flux lies a factor of 2 above the extrapolation of the 5-8 GHz 
flux, and the NVSS image suggests spatial extension on ar- 
cmin scales (a factor of ~ 2 larger than seen at 5-8 GHz). 
This suggests that NVSS is detecting low frequency emission 
from plasma that is emitting a steep radio spectrum, possibly 
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FIG. 2. — Hardness ratio \i2 adaptively binned to reach a signal-to-noise of 5 (after background subtraction) with contours from the 5 — 6<r adaptively smoothed 



X-ray image (contour start from c\ = 10 s ergs cm 



' pixel 1 , increasing by a factor of 1.4 between each contour). 



indicating that it has suffered synchrotron aging. A reasonable 
upper limit on the current radio power can be derived by tak- 
ing the 1 .4 GHz luminosity, and using the arguments of B ick- 
nell et al. (1998) to convert it into a kinetic luminosity. Tak- 
ing the smallest reasonable value of their conversion parameter, 
K1.4 > 10~ 12 , we estimate an upper limit on the current kinetic 
power of Lkin < 7 x 10 43 ergs -1 . 

Comparing the time-averaged source power (derived the X- 
ray cavities) to that derived from the radio luminosity (which is 
equivalent to the source power averaged over the synchrotron 
cooling time of the 1.4 GHz electrons, t coo i <J 10 7 yrs), one 
infers that either this source has faded in kinetic luminosity by 
an order of magnitude or more, or that the magnetic field in 
the lobes is considerably out of equipartition. Since the thermal 
pressure is close to the equipartition pressure estimated by T94, 
we favor the first possibility. Given the uncertainties in these 
arguments (especially in «i 4, for which we chose a conserva- 
tive value), the source could easily have faded by more than 
an order of magnitude. Indeed, the fact that the average power 
is in the realm of FR-II radio galaxies, while morphology and 
current radio luminosity qualify it as an FR-I, leads us to specu- 
late that PKS2354-35 is an example of an FR-II source that has 
faded into an FR-I source on a timescale of less than 10 8 yr. 

The apparent offset between the cluster center and the cen- 
ter of the cavities and the asymmetric brightness distribution 
through the equatorial regions of the hour-glass structure may 
be evidence for bulk ICM motions. In particular, the morphol- 
ogy suggests a bulk flow in a NE direction which might further 
squeeze the outward moving compression wave from the radio 
galaxy, and sweep back the cavity structure. We note that the 
SW ridge is rather cool and thus cannot be a strong shock re- 
sulting from the interaction of a bulk flow with the radio galaxy. 
Hydrodynamic simulations are required to investigate this sys- 
tem further. 



4. CONCLUSIONS 

We have presented CXO observations of Abell 4059. While 
the ICM appears smooth and relaxed on large scales, it shows 
complex morphology in the core region which is likely the re- 
sult of interaction between the ICM and the central FR-I radio 
galaxy PKS2354-35. As was suggested by HS, PKS 2354-35 
appears to have inflated two large cavities in the ICM. Together 
with a central bar-like structure, these cavities produce an hour- 
glass like morphology which can be readily understood as being 
due to a radio cocoon expanding into the ICM. While clear cor- 
respondence exists between the NW cavity and the NW radio 
lobe, the SE cavity is much larger than the SE lobe, suggesting 
that this could be a 'missing link' between cavities with and 
without visible radio lobes. 

The absence of sharp edges in the brightness images and of 
large temperature jumps implies that PKS2354-35 is not driv- 
ing a strong shock into the ICM. We suggest that it is in the 
weak-shock/compression-wave phase identified in the hydrody- 
namic simulations of RHB. Dynamical estimates give a time av- 
eraged kinetic source power of at least Lkm J> 3 x 10 44 erg s -1 , 
while estimates based on the current radio luminosity indicate a 
source power of L\^ u ^ 7 x 10 43 ergs -1 . We suggest that this 
source has faded by a significant amount (and possibly from an 
FR-II phase) during the past 10 s yr. 
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